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Abstract 



Nuclear transparency in the (e, e'p) reaction for 135 <Tp< 800 MeV is inves- 
tigated using the distorted wave approximation. Calculations using density- 
dependent effective interactions, both empirical and theoretical, are compared 
with phenomenological optical potentials. We find that nuclear transparency 
is well correlated with proton absorption and neutron total cross sections and 
that calculations using density-dependent effective interactions provide the 
best agreement with data. Nuclear transparency calculations are compared 
with recent electron scattering data for < 2 (GeV/c)^. For Tp < 200 MeV 
we find that there is considerable sensitivity to the choice of optical model and 
that the empirical effective interaction provides the best agreement with the 
data, but remains 5-10% low. For Tp ^ 300 MeV we find that there is much 
less difference between these models, but that the calculations significantly 
underpredict transparency data and that the discrepancy increases with A. 
The differences between Glauber and optical model calculations are related to 
their respective definitions of the semi-inclusive cross section. By using a more 
inclusive summation over final states the Glauber model emphasizes nucleon- 
nucleon inelasticity, whereas with a more restrictive summation the optical 
model emphasizes nucleon-nucleus inelasticity; experimental definitions of the 
semi-inclusive cross section lie between these extremes. 
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I. INTRODUCTION 



The transparency of nuclear matter to the propagation of an intermediate-energy nu- 
cleon is of fundamental importance to the interpretation of many nuclear reactions. The 
elastic scattering of intermediate-energy nucleons, here defined by the range 100 ^ Tp < 800 
MeV, is usually described using a complex optical potential whose real part produces refrac- 
tion and whose imaginary part produces absorption of the flux contained within the elastic 
channel. That flux reappears, of course, in other channels which include inelastic scattering, 
knockout, and other reactions. Phenomenological optical potentials are obtained by fitting 
the parameters of some hypothetical function to elastic scattering data; sometimes absorp- 
tion or total cross section data are included in the analysis also. Excellent fits to the data for 
a particular target and energy can usually be achieved, but the fitted parameters often vary 
erratically with respect to either mass or energy. Global analyses of more extensive data sets 
which impose smooth dependencies on energy and/or mass, usually at the expense of local 
fit quality, are expected to produce more realistic optical potentials. The optical potential 
is then presumed to represent the nucleon wave function within the medium with sufficient 
accuracy for the analysis of other reactions involving one or more nucleons of similar energy 
in either the initial or the final channel, such as {e,e'N) or {p,p'N). 

However, since elastic scattering is determined by asymptotic phase shifts, any potential 
which produces the same set of phase shifts will predict the same elastic scattering even if the 
interior wave functions are quite different. Thus, various phase-equivalent optical potentials 
can produce significantly different predictions when employed to analyze other reactions. 
Alternatively, we have developed an empirical effective interaction (EEI) [0-^] based upon 
the density dependence predicted by nuclear matter theories of the effective interaction, but 
where the parameters are fitted to inelastic scattering data. Transition densities fitted to 
electroexcitation data are used to minimize uncertainties due to nuclear structure. Both 
elastic and inelastic potentials are obtained by folding the same density-dependent effective 
interaction with transition densities using the local density approximation. Since both the 
distorted waves and the transition potentials depend upon the same interaction, the fitting 
procedure involves a self-consistency cycle which converges quite quickly. Several inelastic 
transitions, possibly among several targets, are fitted simultaneously, producing an empirical 
effective interaction (EEI) that depends upon local density but which is independent of 
target. Elastic scattering data can be analyzed also, but the interaction fitted to inelastic 
scattering data usually produces good elastic predictions whether or not elastic data are 
included. 

Basing the phenomenology of the effective interaction primarily upon inelastic scattering 
offers several advantages over the determination of the optical potential via elastic scatter- 
ing. First, the various radial shapes of transition densities provide differential sensitivity 
to the density dependence of the effective interaction, whereas the elastic optical potential 
requires only a global average over density. Second, since the distorted waves are deter- 
mined self-consistently using optical potentials constructed from the same interaction that 
drives the inelastic transitions, the inelastic observables depend upon overlap integrals and 
are sensitive to the wave functions in the nuclear interior, thereby helping to determine 
the interior optical potential. Elastic scattering, on the other hand, depends on asymptotic 
phase shifts and all potentials which share the same asymptotic wave functions predict the 
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same elastic scattering despite differences in the interior. Therefore, provided that a consis- 
tent description of both elastic and inelastic scattering emerges from the phenomenological 
analysis of the effective interaction, the resultant optical potential should represent a more 
realistic description of the nuclear interior than models restricted to simple geometries, even 
if the latter yield better fits to the elastic scattering data alone. Finally, the EEI model 
requires far fewer parameters than traditional global optical models. 

Single-nucleon knockout by electron scattering, (e, e'A^) reactions, provide important 
tests for models of nucleon propagation. Since the nucleus is practically transparent with 
respect to the electron beam, the ejectile can originate from anywhere within the nuclear 
volume. Similarly, the smaller number of strongly absorbed wave functions gives (e, e'N) 
reactions better sensitivity to the interior than {p,N) or {p,p'N) reactions. Hence, {e,e'p) 
reactions view the optical potential quite differently from {p,p') reactions and there is no 
guarantee that the phenomenological potential that provides the best fit to proton elastic 
scattering data will also provide the best description of (e, e'p) data. 

Quantitative analysis of missing-momentum distributions or spectroscopic factors for ex- 
clusive A{e, e'p)B reactions to discrete states of the residual nucleus depends upon accurate 
knowledge of absorption and distortion by final-state interactions (FSI) 10-0. Most ex- 
periments of that type have been performed using relatively low ejectile energies, typically 
Tp ^ 135 MeV, where the EEI method is less successful. More recent experiments at Mainz, 
with Tp 200 MeV, and future experiments at CEBAF should be less sensitive to FSI 
uncertainties. Nevertheless, it is important to test those models of nuclear transparency 
directly. 

The first measurements of nuclear transparency using electron scattering were made at 
MIT by Garino et al. Il^jll] for Tp 180 MeV by comparing inclusive cross sections with 
measurements of single-nucleon knockout cross sections that sampled the reaction cone and 
which were integrated with respect to missing energy. These semi-inclusive (e, e'p) data 
were analyzed using a correlated Glauber approximation by Pandharipande and Pieper 



12|, who found that Pauli blocking and short-range correlations are required to obtain 



transparencies large enough to reproduce the data. However, the Glauber model [T^ is a 



high-energy approximation which postulates linear trajectories and hence is not expected 
to be particularly accurate for Tp ^ 180 MeV. For proton elastic scattering the Glauber 
model is generally considered adequate for Tp ^ 800 MeV [|1^,0. Although integrated 



quantities, such as nuclear transparency, are probably less sensitive to the details of final- 
state interactions such that the eikonal approximation may be sufficiently accurate at lower 
ejectile energies, the lower limit of the Glauber approximation to (e, e'p) has not yet been 
established. 

At high Q^, Brodsky [|1^ and Mueller [|1^ have predicted that nuclear transparency 
might be significantly enhanced by the phenomenon of color transparency (CT), in which 
a Fock component of the nucleon wave function that is smaller in configuration space than 
the complete nucleon is ejected as an effectively color-neutral object that propagates with 
reduced interactions and increased attenuation length. Early calculations based upon the 
Glauber approximation with energy-independent nucleon-nucleon cross sections and various 
models of the hadron formation length for point-like configurations predicted substantial 
enhancement of nuclear transparency for the A{e,e'p) reaction at high [[18|-[21[]. Nuclear 
transparency data using the A{e,e'p) reaction for 1.0 < < 6.8 (GeV/c)^ have recently 
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been obtained at SLAC by the NE18 collaboration p2| , p3| , but no definitive signal for CT was 
discerned |24]. On the other hand, the data were limited by statistics and new data of much 



better accuracy are expected from CEBAF soon, although the maximum will be smaller. 
However, more refined calculations, including the energy-dependence of the nucleon-nucleon 
interaction |T^, finite range effects [^, Fermi motion [^,^], and coherency constraints 
||28|| , suggest that the onset of color transparency is rather slow and cannot be seen in these 
types of experiments unless ^ 30 (GeV/c)^. Hence, the CT signal in semi-inclusive 
y4(e, e'p) cross sections is expected to be small at SLAC kinematics and even smaller at 
CEBAF kinematics. It has been suggested that for discrete states of the residual nucleus 
the attenuation coefficient for single nucleon knockout [|T^ or the asymmetry between missing 
momenta parallel versus antiparallel to the momentum transfer might be more sensitive 
to color transparency at intermediate energies, but recent studies of the accuracy of the 
Glauber approximation suggest that ^ 2 (GeV/c)^ would be required to employ that 
model 



In this paper we report calculations of nuclear transparency for 135 < Tp < 800 MeV 
using a distorted wave approximation (DWA). We compare calculations based upon the EEI 
model and a relativistic effective interaction, known as IA2 [0,0 , with a global optical po- 
tential from Dirac phenomenology [Q. For low energies we also consider several traditional 
nonrelativistic optical models. The calculations are compared with the MIT data 1 10,11] 
and with the low-energy data from SLAC experiment NE18 p2| , ^ . The model is presented 
in Sec. |I|. Predictions of proton absorption and neutron total cross sections using these 
optical potentials are compared with the available data for these closely related quantities 
in Sec. |T|. The DWA results for nuclear transparency are presented in Sec. In Sec. 
we compare our results with those obtained using the Glauber model. We find that the 
difference between these approaches arises from different definitions of the semi-inclusive 
cross section, with Glauber model being more inclusive and the optical model less inclusive 
than the experimental semi-inclusive cross section for intermediate-energy ejectiles. Finally, 
our conclusions are summarized in Sec. [Vl . 



II. MODEL 

A. Definition of Nuclear Transparency 

The distorted spectral function is obtained experimentally by dividing the differential 
cross section by the off-shell electron-proton cross section, cXep, according to the ansatz 

^^^^^^ = ^-^vS''{E^.V^^^). (1) 

where is a kinematical factor. Final-state interactions between the ejectile and the resid- 
ual nucleus make the distorted spectral function S^{Era-,Vm-,v'^ depend upon the ejectile 
momentum p', and on the angle between the initial and final nucleon momenta, whereas 
the (undistorted) spectral function would depend only on Em and pm- Thus, the distorted 
spectral function depends upon the kinematical conditions and is different for parallel and 
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perpendicular kinematics, for example. Furthermore, the dependence of upon the elec- 
tron energy that arises from Coulomb distortion and from the properties of the electron 
current has been left implicit. 

Due to the distribution of initial momenta, the quasifree single-nucleon knockout strength 
is spread over a Fermi cone whose opening angle is approximately Op — tan~^ kp/q and over 
a range of missing energies that includes the binding energy for the deepest orbital. However, 
final-state interactions broaden these distributions and transfer some of the flux into more 
complicated final states. Hence, nuclear transparency can be loosely defined as the ratio 
between the coincident (e, e'p) cross section and the inclusive quasifree electron scattering 
cross section, where the concidence cross section is integrated over the Fermi cone and over 
the range of missing energy populated by direct knockout. 

Optical models describe the loss of flux from the initial single-nucleon channels as ab- 
sorption. Hence nuclear transparency is calculated in the optical model as the ratio 

JdE^Sd^PmW{E^,Pm,p')SP^{Em,Pm) ^' 

where in the numerator the distorted spectral distribution includes final-state interactions, 
and depends upon the ejectile momentum, whereas in the denominator the undistorted 
spectral distribution depends only upon the missing energy and momentum. The weight 
factor, w{EmiPm.-,p')^ represents the experimental acceptance and distinguishes, for exam- 
ple, between parallel or quasiperpendicular kinematics; hence we distinguish various types of 
transparency functions using a subscript w signifying the appropriate acceptance function. 
Clearly distortion and transparency depend upon the ejectile energy. However, since the 
distorted spectral function must be evaluated for the appropriate kinematical conditions, its 
dependence upon electron-scattering kinematics remains implicit in this definition. There- 
fore, to complete the definition of transparency we would need to specify the kinematical 
conditions of interest more completely, using the same integration regions and weighting fac- 
tors for both numerator and denominator, where the appropriate weighting factors depend 
upon the kinematics of the experiment. 

The simplest situation arises when (a;, q) are held constant, for which the nuclear trans- 
parency 71 may be defined as 

_ J dEm J d9 sine S^{E„„Pm,p\9) 
^~ JdE^JdesmeSP^{Err^,Pm) 

where 9 is the angle between the ejectile momentum p' and the momentum transfer q. 
Furthermore, it is simplest to require that the ejectile kinetic energy be constant in the 
barycentric frame so that final-state interactions can be evaluated for a unique ejectile en- 
ergy. Also note that the weight factor should be expressed in terms of the center-of-mass 
system, where the ejectile momentum is constant, but for large q the laboratory momentum 
changes little over the Fermi cone so that the distinction between the lab and cm angles 
matters little. However, for experiments using a narrow acceptance in ou with a large ac- 
ceptance in ejectile kinetic energy Tp, it is necessary to employ energy-dependent optical 
potentials. Alternatively, for parallel kinematics we define 
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/ dE^ J dp 

" ~ JdE^Jdprr^pl SP^{E^,p^) ^ ^ 

but must recognize that the ejectile momentum p' = pm + 5' is correlated with p^ and that 
the electron scattering kinematics again remain implicit. In this case it is also clearest to 
require the invariant mass for each E„i to be constant despite the concomitant variation of 
electron kinematics. 

We have verified that when the optical potentials are nullified, the distorted wave calcula- 
tions result in unit transparency. However, it is important to recognize that these definitions 
of nuclear transparency can produce deviations from unity, with either sign, even when the 
optical potential is purely real because weighting functions based upon the plane-wave im- 
pulse approximation do not account for refractive effects. Although one might be tempted 
to divide out the refractive effect of the real potential by using a modified definition 

^ _ / dEjn J d^Pm w{Em, Pm, p')S^{Em, Pm, PQ , . 

J dEm J d^Pm w{Em, Pm, p')S^{Em, Pm, p') 

where S^{Em, Pm, p') is the distorted spectral function for the real part of optical poten- 
tial, our original definition, Eq. (^), conforms more closely to the customary experimental 
definition. The difference between these approaches is generally greatest when the summa- 
tion over missing energy is restricted to a single state, or subshell. For quasiperpendicular 
kinematics, the refractive effects are minimized when measurements are made on both sides 
of q for closed shell targets so that the opposing effects of spin-orbit distortion are ap- 
proximately balanced for both sides of q and for both spin-orbit partners. For parallel 
kinematics a net attractive (repulsive) real central potential shifts the missing momentum 
distributions to smaller (larger) pm and enhances (reduces) the peak values of S^{Em,Pm) 
for positive (negative) missing momenta, but the net effect on transparency tends to bal- 
ance when the integration over p^ is symmetric. Since the shifts in the peak positions are 
also sensitive to spin-orbit distortion, the attenuation factors for spin-orbit partners can be 
different. Nevertheless, numerical studies show that for calculations using symmetric ranges 
of missing momenta the ratio T^/'T^ remains within a few percent of unity for both parallel 
and quasiperpendicular kinematics for A>12 and 135 <Tp < 800 MeV. Alternatively, one 
can assess the effects of refraction by comparing transparency calculations with and without 
the real parts of the optical potential. Thus, over the ranges of A and Tp considered, we find 
that \5T\/T < 7% when the real parts of the optical potential are eliminated and conclude 
that the present definition does in fact provide an unambiguous measure of transparency 
that is rather insensitive to refraction. Therefore, we prefer to employ Eq. (^, despite its 
slight ambiguity between refraction and absorption, because it is closer to the experimental 
quantity of interest and the refractive effects do not significantly affect its interpretation in 
terms of absorption. Furthermore, we also find that 



77, -T, 



r, 



< 6% (6) 



for A > 12 and 135 < Tp < 800 MeV; hence, we expect our calculations to be quite 
insensitive to the small range of missing momentum parallel to q that is inevitably accepted 



6 



when performing measurements of 71. Therefore, we compare calculations of 71 directly to 
the experimental data without attempting to simulate the complete experimental acceptance 
functions. 



B. Spectral Function 

In the independent-particle model (IPM) the spectral functions take the forms 
Sipm{E 

) = H Saga{Pm)S{Em - E^) (7) 
a 

Sfp^iEm, Pm, P) = E SaQ^^iPm, p')S{Em - E^) (8) 



where E^ and 5"^ are the missing energy and spectroscopic factor for orbital a, 

2 



Qa{,Pn 



(9) 



is the momentum distribution arising from a single-particle wave function 0Q,(r) normalized 
such that 

An J dpm pin QaiPm) = 1 , (10) 

and 0aiPm,p') is the corresponding distorted momentum distribution. Hence the nuclear 
transparency for quasiperpendicular kinematics reduces to 

^ ^ Eg Sg J de sin e (p^, p', 9) 
^ HgSgj d9 sin 9 (pm) 

Although interactions fragment the low-lying valence hole strength and populate the 
continuum for large missing energy and momentum, nuclear transparency is much more 
sensitive to final state interactions than to details of the spectral distribution. For a given 
orbital, attenuation depends more strongly on ejectile energy than on missing energy, so 
that the spreading of the hole strength is not expected to appreciably affect its contribution 
to the integrated yield, especially for experiments which select a relatively narrow range of 
ejectile energy. Similarly, although variations of the bound-state wave function which change 
the rms radius affect the missing momentum distribution, such variations have little effect 
on the integrated yield and tend to cancel in the ratio used for transparency. Therefore, 
for transparency calculations it is sufficient to employ IPM spectral functions, although 
correlations must be included to describe the distorted spectral distribution in detail. 

For A < 16 we used Woods-Saxon single-particle wave functions, whereas for heavier nu- 
clei bound-state wave functions were obtained using the Skyrme-Hartree-Fock model based 



upon the interaction designated Z^r The single-particle energy spectra for each tar- 

get were shifted to obtain the correct separation energies. For each orbital the final-state 
interactions were evaluated using ejectile energies based upon u and the shifted Hartree- 
Fock separation energies. Note that the results are rather insensitive to the details of the 
bound-state wave functions, but are quite sensitive to the choice of optical model. 
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C. Distorted Wave Approximation 



The distorted wave approximation (DWA) for the electromagnetic transition amphtude 
that governs the single-nucleon knockout reaction A[e, e'N)B can be expressed in the form 
31 



^ = (12) 

where the electron and nuclear currents are 

= J d'r e-'^'-^rf{rhM^) (13a) 
J-j^(q') = I d'r e^'i'-V^f (r)r,^f (r) (13b) 

and where is the vertex operator for the nucleon current. In these expressions the electron 
wave functions relative to the target of mass rriA are denoted by the spinors ipf and tpj for 
the initial and final states, respectively. At this stage we leave implicit the dependence of the 
nuclear current upon the ejectile kinematics and the state of the residual nucleus. Since it 
is more convenient to express the ejectile wave functions relative to the residual nucleus 



of mass niB, the radial scale is adjusted by means of the reduced momentum transfer ||37 
q' = c(mB/mA- 

If we assume that a virtual photon with momentum q' is absorbed by a single nucleon 
with initial momentum p, the nuclear current at position r becomes 



where the single-nucleon wave function is the amplitude for removing a nucleon from the 
initial state of target A and reaching the final state of residual nucleus B, such that 

0(p) = (S|a(p)|A). (15) 

The distorted wave X^~''*(p') p") is the amplitude that the ejectile with initial momentum 
p" = p + q' emerges from the nuclear field with final momentum p'. In coordinate space 
these wave functions are expressed as 

0(r) = / ^ e'">(p) (16a) 
X(p',r)= /|^e^""-W,P")- (16b) 
Thus, the nuclear current becomes 

:7-j"(p', q) = / x^"^*(p , p + q')r,(p + q', p)0(p) (17) 

where q' is the local momentum transfer supplied by the electron. 
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Since nuclear transparency depends upon the nuclear final-state interactions and not the 
electronuclear initial-state interactions, distortion of the electron wave function should either 
be included in both numerator and denominator of Eq. (Q) or excluded from both. To a good 
approximation, Coulomb distortion can be described as a shift of the effective momentum 
transfer and a focussing factor which increases the virtual-photon flux p8|j39|] . Since these 
aspects of Coulomb distortion have similar effects upon both exclusive and inclusive electron 
scattering, it is reasonable to omit Coulomb distortion for both. In the absence of Coulomb 
distortion, the electron current is proportional to 5'^(q' — q), so that the nuclear current 
can be evaluated for a unique value of the momentum transfer obtained from asymptotic 
kinematics. Therefore, we obtain 



^j"(p',q) 



d'r e^^V-^*(p',r)r^(p',p'-q)0(r) 



(18) 



where the vertex function has now been reduced to a matrix, acting on nucleon spins, whose 
elements are evaluated using effective kinematics. 



We used the ccl vertex function of de Forest with nucleon form factors from model 3 

Current conservation was enforced at the one-body level 



of Gari and Kriimpelmann pl| , |42 



by eliminating the longitudinal in favor of the charge component. However, it is important 
to note that the transparency calculations are quite insensitive to these choices for the vertex 
function and, in fact, are quite insensitive to the electron-scattering kinematics also. The 
most important variables are the ejectile energy and the choice of optical model. 



D. Optical Models 

Distorted waves were obtained from solutions to a Schrodinger equation of the form 

(V' + kl- 2fi^U^)Uka, r) = (19) 
where relativistic kinematics are incorporated by interpreting ka as the exact relativistic 



wave number and fia as the reduced energy for channel a |^3[. We assume that the optical 
potential can be reduced to local form and that nonspherical components may be neglected. 
Thus, the optical potential takes the form 

Uair) = f/f (r) + f/f (r) + f/„^^(r)L ■ cr (20) 

1 dF^^ 

where is the Coulomb potential, U^' is the central potential, and U^''^ = — is 

r or 

the spin-orbit potential for exit channel a. Although in principle the optical potential for 
each exit channel depends upon the structure of the residual nucleus, we employ a mean 
field approximation for single-nucleon knockout in which there is only a small kinematic 
dependence upon the ejectile energy. To accomodate models which include nonlocality 
corrections, we identify the distorted wave function as 

Xa{ka,r)=P{r)Uka,r) (21) 

where P{r) is a Percy factor which is unity for local models, approaches unity at large 
distances for nonlocal models, and which may be complex. 



9 



1. Nonrelativistic Woods-Saxon Potentials 



The potential fitted by Schwandt et al. to cross section and analyzing power for 
proton elastic scattering for A > 40 and 80 < Tp < 180 MeV is commonly employed for 
knockout analyses. It may be used for either proton or neutron scattering since it contains a 
parametrization of the symmetry potential. Although the Schwandt potential does not bear 
extrapolation in either mass or energy well, it is nevertheless often used for lighter nuclei, 
sometimes even for mass-12. A potential developed by Abdul- Jalil and Jackson P5|,^ for 



A ^ 12 and 50 < Tp < 160 MeV has sometimes been used for knockout studies, but in our 
opinion its description of proton scattering data is unsatisfactory. Alternatively, the p+^^C 



potential of Comfort and Karp for Tp < 185 MeV is preferred. Unfortunately, a global 
nonrelativistic optical potential for A > 12 and a broad range of energy does not appear to 
exist. For that we must appeal to Dirac phenomenology (DP). 

Many analyses of single-nucleon knockout also include a Perey nonlocality correction of 
the form mm 



P(r) 



1-1/2 

(22) 



where the central potential is separated into real and imaginary parts denoted U^{r) = 
V'-'{r) + iW'"{r). The nonlocality parameter is typically chosen as /5 = 0.85 fm, based upon 
the original analysis of neutron scattering for T„ < 20 MeV. Although the applicability of this 
simple prescription has not been established for Tp > 100 MeV, when using nonrelativistic 
Woods-Saxon potentials, such as the Schwandt or the Comfort and Karp models, we conform 
to standard practice by including the Perey factor. Since the Perey factor is less than unity in 
the interior and equal to unity outside the potential, its effect is to reduce the transparency, 
especially for interior orbitals. 



2. Dirac Phenomenology 



Suppose that a four-component Dirac spinor 



^+(r) 



where and '0_ are two-component Pauli spinors for positive and negative energy compo- 
nents, satisfies a Dirac equation of the form 



[a ■ p + P{m + S)]"^ = {E -V - 1/^)^ 



(23) 



with scalar and vector potentials S and V. Upon elimination of the lower components, an 
equivalent Schrodinger equation of the form 



+ A;2 - 2^ (U^ + + U^^L 







(24) 



can be obtained, where is related to by the Darwin transformation 
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i)+ = B^/^(j) (25a) 
S -V -V^ 

B = l + -—- . (25b) 

h + m 

The Schrodinger solutions are phase-equivalent to the Dirac solutions in the sense that the 
asymptotic phase shifts, and hence observables for elastic scattering, are the same. However, 
the Dirac wave function is modified in the interior by a nonlocality factor similar in form to 
the Perey-Buck nonlocality factor, except that it depends upon S—V, which is closely related 
to the spin-orbit potential, rather than upon the central potential. Thus, B can be deduced 
directly from the spin-orbit potential |^|. When used in nonrelativistic calculations, the 



positive energy spinor ?/^+, including the Darwin factor, is identified with the distorted wave 
X- 

Hama et al. |^ produced global Dirac optical potentials for A > 40 and 65 < Tp < 1040 
MeV. The global Dirac optical potential was then extended by Cooper et al. to the ranges 
A > 12 and 20 < Tp < 1040 MeV. Of the several essentially equivalent variations of the 
global potential that were provided, we have chosen the version labelled EDADl. Although 
the scalar -l-vector (SV) model of Dirac phenomenolgy is not unique, and relatively simple 
but arbitrary shapes are employed for the potentials, this work represents the most extensive 
and systematic analyses of proton optical potentials available. The available proton-nucleus 
elastic scattering data are described very well by potentials whose properties vary smoothly 
with both mass and energy. 



E. Local Density Approximation 



Optical potentials can also be obtained by folding the nucleon-nucleon effective inter- 
action with the nuclear density distribution. In recent years it has become clear that the 
intermediate-energy nucleon-nucleon effective interaction depends strongly upon the density 
in the interaction region. Several calculations of the effective interaction in nuclear matter 
have been made following the seminal work of Hiifner and Mahaux |]52|. Jeukenne, Lejeune, 
and Mahaux [p^-|56[| computed the self-energy and the optical potential for < Tp < 160 
MeV using the Reid soft-core potential Brieva, Rook, and von Geramb [^^ -p^ devel- 
oped a Brueckner-Hartree-Fock (BHF) approach and used the Hamada- Johnston potential 
pB| to calculate the pair wave function in nuclear matter, from which a local pseudopo- 
tential was constructed for Tp < 180 MeV using a generalization of the Siemens averaging 
procedure ||6^. Similar calculations using the Hamada- Johnston potential have also been 
performed by Yamaguchi, Nagata, and Michiyama (YNM) who parametrized their 

results for Tp < 200 MeV in Gaussian rather than Yukawa form. The BHF approach was 
refined by von Geramb and collaborators p7|j68|] , who constructed an effective interaction 
based upon the Paris potential |^9|], designated Paris-Hamburg (PH), that is applicable for 
100 < Tp < 400 MeV. Nakayama and Love used the Bonn potential to calculate 
a local pseudopotential that reproduces on-shell matrix elements of the G-matrix. These 
theories are all based upon the Bethe-Goldstone equation and include Pauli blocking and 
self-energy corrections self-consistently. A closely related calculation by Ray [|72[ used a 



coupled channels nucleon-isobar model and Watson multiple scattering theory to calculate 
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a density- dependent t-matrix applicable to nucleon energies above 200 MeV that includes 
some of the effects of pion production. Finally, Furnstahl, Wallace, and Kelly [31,32| have 
developed an effective interaction in a form similar to the EEI model based upon the rel- 
ativistic IA2 model. Density dependence arises from the distortion of Dirac spinors in the 
nuclear medium, which primarily affects the real central interaction, and Pauli blocking, 
which damps the absorptive potential. The model naturally provides stronger density de- 
pendence for inelastic scattering than for elastic scattering, which is needed to describe the 
data. 

All of these calculations predict strong density dependence of the nucleon-nucleon ef- 
fective interaction. The dominant effect for the real central interaction is equivalent to a 
short-ranged repulsive interaction that is proportional to density and nearly independent of 
energy. In the BHF approaches this short-range repulsive interaction arises both from the 
anticorrelation between identical nucleons in nuclear matter and from dispersive effects in 
the self-consistent mean field, whereas in the IA2 approach it arises from spinor distortion in 
the strong scalar and vector mean fields. However, both approaches predict that the density 
dependence of the real central interaction depends slowly upon energy and remains quite 
strong even at 800 MeV ||43|j3^ . The dominant effect for the imaginary central interaction 
arises from Pauli blocking and for both the BHF and the IA2 models gives results similar to 
the familiar Clementel- Villi |]T3[ damping of the absorptive potential, in which the damping 
factor is inversely proportional to the proton energy. Hence this effect is most important for 
low energies. 

Although the qualitative features of the medium modifications are essentially the same, 
the quantitative differences among the various theories are surprisingly large, much larger 
than would be expected from the variations among the underlying nucleon-nucleon poten- 
tials that are employed, suggesting that the approximations required to evaluate the effective 
interaction are not yet under good control. The effect of these differences upon elastic and 
inelastic scattering calculations has been surveyed in a series of papers by Kelly and col- 
laborators in which transition densities measured by electroexcitation are used to minimize 
uncertainties due to nuclear structure and to isolate the effective interaction for detailed 
examination Transition densities which are strong in the interior provide informa- 

tion about the high-density properties of the effective interaction, whereas surface-peaked 
transition densities reveal the low-density properties. The systematic comparison of such 
cases demonstrates quite clearly that the effective interaction depends upon local density 
and that estimates based upon nuclear matter theory have qualitatively correct characteris- 
tics, but that none of the theories presently available is sufficiently accurate for quantitative 
applications to nuclear structure. Therefore, an empirical model of the effective interaction 
was developed in which medium modifications similar to those predicted by nuclear matter 
theory are parametrized in a form suitable to phenomenological analysis of data. The pa- 
rameters are adjusted to reproduce inelastic scattering data for several states in one or more 
targets simultaneously. Empirical effective interactions have been extracted from data for 
several energies in the range 100 <Tp< 650 MeV. For each energy we find that a unique ef- 
fective interaction describes data for several inelastic transitions in a single nucleus and that 
the fitted interaction is essentially independent of target. These findings tend to confirm 
the basic hypothesis of the local density approximation (LDA), namely that the interaction 
depends primarily upon local density and is independent of the detailed structure of any 
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particular target or transition. The fitted parameters also exhibit a relatively smooth energy 
dependence. 

Although slightly better fits to some of the data sets may be found in the original analyses, 
for the present purposes we choose to employ the empirical effective interactions (EEI) 
tabulated by Kelly and Wallace [^]. For this set of interactions, the medium modifications 
are applied to the Franey-Love (FL) |j74| parametrization of the free t-matrix, which is 
available for all relevant energies, and common fitting strategies and constraints were used 
to help smooth the energy dependence of the fitted parameters. These choices are made 
primarily for aesthetic reasons and have very little effect upon knockout calculations. 

In the local density approximation (LDA), the central and spin-orbit potentials become 

U'\r) = - fdq q^Uqr)t^{q,p)p{q) (26a) 

F'-'ir) = - [dq q'Mqr)t"^' {q, p)p{q) (26b) 

where 

p{q) = dr r'^jQ{qr)p{r) (27) 



is the Fourier transform of the ground-state density p. Note that a sum over nucleon or 
isospin indices has been left implicit. To minimize uncertainties due to the nuclear density, 
the proton density was obtained by unfolding the proton charge form factor from the charge 



densities measured by electron scattering and tabulated in Refs. For relatively small 

momentum transfers, charge symmetry ensures that the neutron and proton densities for 
mirror nuclei are very nearly proportional to each other and since the high q properties 
of distorting potentials have very little effect upon knockout calculations, especially for 
integrated strengths, we use p„(r) = Npp{r)/Z. Furthermore, we evaluate the local density 
at the site of the projectile. 

Since the empirical effective interaction was fitted to data for elastic and inelastic scat- 
tering self-consistently without explicit use of a nonlocality correction, no Percy factor is 
used with the EEI model for knockout. Similarly, the IA2 interaction is local by construc- 
tion and does not require a Percy factor either. Some other versions of the LDA do involve 
nonlocality corrections based upon the exchange contribution or upon the momentum de- 
pendence of the effective mass, but are not employed here because none of those models 
provide adequate descriptions of the proton scattering data. 

The potentials which emerge from all microscopic models exhibit much more compli- 
cated radial shapes than posited by the Woods-Saxon model of the optical potential. The 
detailed shapes depend upon the density dependence and range of the effective interaction 
and upon the nuclear density, which especially for light targets is not well approximated 
by the Fermi shape. For energies between about 100 and 300 MeV, for example, the real 
central potential exhibits a characteristic "wine-bottle" shape. Similar shapes also arise 
from the nonrelativistic reduction of either Dirac phenomenolgy or the relativistic impulse 
approximation. Although good fits to elastic scattering data may be achieved with sim- 
plistic models of the potential, artificially simple geometries cannot be justified on more 
fundamental grounds. Furthermore, the missing momentum distributions for discrete states 
do show some sensitivity to the shape of the real central potential. 
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III. COMPARISON OF OPTICAL MODELS 



Integrated quantities, such as nuclear transparency, are much more sensitive to attenu- 
ation of the flux than to distortion of the angular distribution by final state interactions. 
Therefore, in this section we compare data for proton absorption and neutron total cross 
sections with calculations based upon a variety of optical models for 100 < Tp < 800 MeV. 
More detailed analyses of proton elastic and inelastic scattering data can be found in Refs. 
|]T|-0, for example. 

Predicted proton absorption cross sections are compared with data in Fig. 0. Unfortu- 
nately, the available proton data are scarce and of uneven quality |TB|. Nevertheless, the 



EEI model provides accurate predictions for these data, although the 200 MeV interaction 
appears to give results which are systematically low compared to the trends for other ener- 
gies. It is important to remember that the EEI model is dominated by inelastic data, but 
gives good fits to elastic data whether or not they are included in the analysis [[77|. Fur- 
thermore, neither absorption nor total cross section data were included in the analysis, but 
are nevertheless predicted accurately. Dirac phenomenology also provides good predictions 
for proton absorption cross sections, but its calculations for both ^^C and ^°Ca appear to be 
slightly too large. Below about 150 MeV the Schwandt model also agrees with the data, but 
its energy dependence appears to be unreasonable and it begins to diverge from the data 
for higher energies. The earlier version of that model due to Nadasen ||7^ gives a better 
description of the energy dependence of the absorption cross section. 

In Fig. ^ predictions for neutron total cross sections are compared with the high-quality 
neutron total cross section data recently obtained at LAMPF fT^. Neutron total cross 
sections computed from the Schwandt potential are substantially larger than the data. For 
self-conjugate targets Dirac phenomenology provides good predictions over broad ranges of 
mass and energy, but, lacking a parametrization of the symmetry potential, the EDADl 



model fails to reproduce (T„ for N > Z pO[- As for the proton absorption cross section, 
EDADl predictions for the ^^C neutron total cross sections appear to be slightly too large 
in the energy range relevant to existing proton knockout data for discrete states. The EEI 
model also provides good predictions for self-conjugate targets and appears to be more 
accurate than EDADl for light nuclei, such as ^^C. The result at 200 MeV appears to be 
slightly too small, as also observed in the proton absorption cross sections. Below about 
below 300 MeV, the EEI calculations for heavier targets with significant neutron excesses 
also appear to be more accurate than EDADl, but tend to be too high at higher energies, 
particularly for ^^'^Pb. For these EEI calculations the symmetry potential is obtained by 
folding the isovector density with the density-independent isovector interaction from the 
FL t-matrix. Although the isovector interaction has not been calibrated to nucleon-nucleus 
scattering data with the same care as the isoscalar interaction, its contribution is small 
enough that residual errors in theoretical models of that term should not be too serious. On 
the other hand, even though the EEI model provides a good fit to proton elastic scattering 
by heavy nuclei, there may still be some inaccuracy for large A because the model was fitted 
to data for A < 40. 

We find that the IA2 interaction provides accurate predictions for elastic and inelastic 
scattering for 500 MeV protons, where the IA2 interaction is most similar to the empirical 
effective interaction, but that the real-central repulsion of the IA2 model is too strong 
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at lower energies. Nevertheless, Fig. ^ shows that very accurate predictions are obtained 
for neutron total cross sections [^. Similarly, calculations using the relativistic impulse 
approximation, with |^ or without density dependence, also reproduce the the neutron 
total cross section data for ^°^Pb. Thus, it will be of interest to compare IA2 predictions 
with measurements of nuclear transparency at higher energies soon to be made at CEBAF. 

In Figs. |l| and ^ we also show calculations based upon the Paris-Hamburg (PH) interac- 
tion. Although the PH optical potential provides good fits to proton elastic scattering data 
|67| , |82[| , we find that its predictions for integrated cross sections are substantially larger than 



the data. Similarly, PH calculations for proton inelastic scattering to states with surface- 
peaked transition densities also tend to produce cross sections that are too large [Q,^,^ . We 
have argued that for finite nuclei nonlocal corrections to the LDA suppress the interaction 
strength in the surface region, producing smaller inelastic cross sections for surface-peaked 
states and smaller integrated cross sections. This effect is included in the empirical effec- 
tive interaction and accounts in part for the improvement of the EEI with respect to the 
PH model for these quantities. Nevertheless, the PH model provides equal or better fits to 
proton elastic scattering. Thus, there is no guarantee that optical models fitted to elastic 
scattering data alone will provide the best predictions for absorption or total cross sections 
or for nuclear transparency. Phenomenological models, such as EEI, which include data 
that are sensitive to the interior wave function, such as proton inelastic scattering, should 
provide more accurate interior potentials and better predictions for nuclear transparency. 
Furthermore, we have also shown that cross sections for exclusive (e, e'p) reactions to dis- 
crete final states correlate well with absorption or total cross sections for nucleon-nucleus 
scattering H]. 



IV. RESULTS 

A. Comparison with data for Tp k, 180 MeV 

The first measurements of nuclear transparency for (e, e'p) with 150 < Tp < 210 MeV 
were made by Garino et al. JTU| , p]T| using a 780 MeV electron beam, uj = 215 ± 20 MeV, and 



q ^ 610 MeV/c. Measurements were made for four in-plane opening angles between about 
0° and 23° with (p = 180°, which sample a slice through the Fermi cone for this q. The ratios 

= dMe,e')/dn^ (^^^ 

integrated over a wide range of missing energy (up to about 100 MeV) were measured for 
^^C, ^^Al, ^^Ni, and ^^^Ta. Several methods were used to relate the measured ratios TZ{0) 
to the nuclear transparency, but the variations were less than 5%, which can be viewed as 
an estimate of the systematic error. From our point of view, the experimental realization of 
Eq. (|11D would be 

_ jdesmon^^je) 
^ ~ jdesinenp^ie) ' ^ ^ 
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where the range of integration is restricted to the large angle side of q. The data for Em < 80 
MeV are shown in Fig. 0. 

Since this experiment was performed with a relatively narrow acceptance in electron 
energy, the calculations were performed with fixed u = 215 MeV and q = 605 MeV/c. The 
proton spectrometer, on the other hand, accepted a wide range of proton energies centered 
about 180 MeV. Hence, for each shell the laboratory proton energy for parallel kinematics 
was computed as Tp = — Em where Em is the Hartree-Fock single-particle energy. The 
invariant mass for each final state was held constant, so that the laboratory proton energies 
for nonparallel kinematics are slightly smaller. The optical potential for each shell was 
calculated using the proton energy for that shell. 

Unfortunately, the EEI interaction is only available for discrete energies and the param- 
eters do not vary as smoothly as one might like because those interactions were obtained 
by fitting uncorrelated data sets independently. No attempt has yet been made to impose 
a smooth energy dependence upon the empirical effective interaction. Nevertheless, the in- 
teractions do vary smoothly enough that calculations made using the 180 MeV or the 200 
MeV interaction are quite similar and the small differences between them do not affect the 
conclusions. Furthermore, a large part of the energy dependence of the effective interaction 
arises from the exchange contributions, which were evaluated using the kinematics of each 
orbital. Therefore, we chose to use the fit to data for 180 MeV protons, representing the 
center of the proton energy distribution. Similar calculations using the 200 MeV EEI inter- 
action, fixed u rather than fixed Tp, and a other few technical differences, were presented in 
Ref. 0. Although slightly better agreement with the data was obtained using that slightly 
less absorptive interaction, the small difference between the two calculations demonstrates 
that the results are relatively insensitive to uncertainties in the EEI parametrization and 
ambiguities in the prescription for acceptance averaging. 

We have studied the sensitivity of nuclear transparency to the optical model by compar- 
ing the EEI model with Dirac phenomenology, version EDADl. Although the data extend in 
both energy and mass beyond their ranges of applicability, we also show calculations based 
upon the Schwandt and the Comfort and Karp (CK) optical potentials for comparison with 
other authors who have used those models. Ejectile wave functions for the Schwandt and 
CK potentials include Perey factors, Eq. (|22|) , with (3 = 0.85 fm, whereas wave functions 
based upon Schrodinger-equivalent potentials from Dirac phenomenology include the Dar- 
win factor, Eq. (|25|) . Nonlocality corrections are not needed for the EEI model because both 
elastic and inelastic scattering are fitted self-consistent ly in that model. These calculations 
are compared with the data in Fig. |^. We find that the EEI model provides a good de- 
scription of nuclear transparency, whereas considerably smaller transparencies are obtained 
with either the Schwandt or the EDADl potentials. These results are consistent with the 
observation in Sec. |TTl| that the EEI model provides more accurate predictions for proton 
absorption and neutron total cross sections also. It is clear that the EDADl model is too 
absorptive at these energies, so that spectroscopic factors using it will be artificially large 
to compensate for excessive attenuation. 

Similar calculations have been performed by Ireland et al. |8^ using the Schwandt po- 
tential and they obtained larger transparencies which agree better with the data than do 
our ostensibly similar calculations with the same potential. However, that analysis suffers 
from several defects. First, the kinematics were artificially altered so that knockout from 
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every orbital was assigned the same ejectile energy, 180 MeV, despite the wide range of en- 
ergies covered by the experiment. Thus, they used the Schwandt potential at 180 MeV even 
though, for the same electron kinematics, protons ejected from less bound orbitals actually 
emerged with energies beyond the range of that optical model. Second, their calculations 
used a scattered electron energy that is 30 MeV higher than the center of the experimental 
acceptance. Third, the DWIA cross section was computed using a current operator based 
upon nonrelativistic reduction to second order in q/rriN-, but the distorted momentum dis- 
tribution was obtained by dividing the cross section by (Ted, the cross section for the ccl 
current operator. The inconsistency between the numerator and the denominator in their 
application of Eq. leads to a spurious enhancement of the transparency above unity for 
a plane-wave calculation. However, none of these defects appears to be sufficient to explain 
the discrepancy between the two calculations. Nevertheless, we consider the apparent agree- 
ment between the data and the Schwandt calculations of Ref. |^ to be a fortuitous result 
of an incorrect calculation. 

Phenomenological optical models, nonrelativistic or relativistic, which are fitted only to 
elastic scattering data may fit that data well but still fail to predict nuclear transparency cor- 
rectly because such analyses are not sensitive to the interior wave function. The EEI model 
is much more sensitive to the interior wave function because in fitting inelastic scattering 
data it requires consistency between distorted waves and inelastic transition amplitudes. 
Therefore, the EEI model provides a more accurate prediction of nuclear transparency and 
should also provide more accurate spectroscopic factors. Clearly it will be of interest to 
obtain transparency data for a wider range of energies and such studies are planned for 
CEBAF. It would also be of interest to obtain comparable data at lower energies where 
NIKHEF has performed its extensive survey of spectroscopic factors for complex nuclei 0, 
but it appears that the lower-energy regime will soon be abandoned. 



B. Comparison with data for Tp 650 MeV/c 

Nuclear transparency data for 1.0 < < 6.8 (GeV/c)^ have recently been obtained at 
SLAG by the NE18 collaboration [p2| , p3| . In Fig. ^ we compare calculations with the data 
for = 1.04 (GeV/c)^ and lj = 0.625 GeV. The angular integrations were symmetric with 
respect to q. Since the dependence of nuclear transparency on proton energy is quite slow 
for these energies (see following section), we performed the calculations using a fixed proton 
energy of 650 MeV for which both the EEI and IA2 interactions are available. The electron 
kinematics were then computed using u = E^+Tp for each shell. At this energy we find that 
the EEI, IA2, and EDADl potentials all give practically the same results — the variation 
due to choice optical potential is much smaller at 650 than at 180 MeV. However, although 
the calculations are fairly close to the data for ^^G, the data vary much less with A than do 
the calculations. O'Neill et al. report that the data can be fitted with a power law of the 
form y4~" with a = 0.18 at = 1.0 (GeV/c)^ and with slightly larger values of a at higher 
Q^. However, the calculations are much closer to the characteristic a = 1/3 behavior that 
would be expected for an eikonal model with constant attenuation length. 

The experimental values of T± reported by the NE18 collaboration [P^,P3] include cor- 



rections for the portion of the model spectral functions which would fall outside their ac- 
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ceptances. These correction factors are 1.11, 1.22, and 1.28 for carbon, iron, and gold, 
respectively. The calculations of Nikolaev et al. appear to be in good agreement with 
the A-dependence of the data for Q"^ ^ 1.0 (GeV/c)^, but their figure apparently employs 
a preliminary analysis of the data or neglects these correction factors resulting in smaller 
experimental values. However, elimination of the correction factors would not be enough 
to bring the data into agreement with our calculations. As discussed in Sec. 0, the most 
important difference between our calculations and those of Nikolaev et al. can be traced to 
differences between the definitions of the semi-inclusive cross section from which transparen- 
cies are calculated. Nevertheless, it would be of interest to compare proton absorption and 
neutron total cross sections computed with their Glauber model with nucleon-nucleus data 
for this energy regime. 



C. Energy Dependence of Nuclear Transparency 

The energy dependence of nuclear transparency is examined in Fig. ^, which com- 
pares calculations using the EEI, IA2, and EDADl potentials for representative light (^^C), 
medium (^^Ni), and heavy (^'^'^Pb) nuclei. The calculations were performed for proton en- 
ergies marked by symbols, where circles are for ^^C, squares are for ^^Ni, and diamonds are 
for ^°^Pb. The calculations were made for constant (tu, q) kinematics using a beam energy 
of 2.1 GeV, although the results are practically independent of the electron energy. The 
energy transfer for each shell was taken to be = Em + Tp, where Tp is the laboratory 
kinetic energy for zero recoil and Em is the Hartree-Fock energy for the shell. The angular 
integrations were symmetric with respect to q. 

The proton energies were selected according to the availability of EEI and/or IA2 inter- 
actions. The small kinks in the EEI calculations between about 180 and 200 MeV probably 
indicate the degree of model dependence that arises when fitting independent data sets 
without imposing a smooth energy dependence. These kinks are clearly correlated with the 
corresponding neutron total cross section calculations shown in Fig. ^ which suggest that 
the 200 MeV empirical effective interaction is not quite sufficiently absorptive. On the other 
hand, the EEI calculations of neutron total cross sections for heavy nuclei tend to be a little 
too high at 180 MeV. Hence, the best estimate of the transparency between 180 and 200 
MeV probably lies between the calculations shown in Fig. ^ with an uncertainty comparable 
to the difference between them. 

For energies below about 200 MeV, the EEI model is more transparent than any of the 
other optical models we have examined and provides the most accurate predictions for the 
MIT data. Nevertheless, those data still remain about 5-10% above the EEI calculations. 
At higher energies the variation among optical potentials tends to decrease, although for 
medium and heavy nuclei the IA2 calculations are more transparent than either the EEI 
or EDADl models. This behavior is also clearly correlated with the neutron total cross 
section calculations shown in Fig. |^, where for medium and heavy nuclei the IA2 interaction 
provides accurate predictions while the EDADl and EEI calculations are similar and are 
both slightly too large. 

A series of ^^C(e, e'p) experiments performed at MIT-Bates in parallel kinematics near 
the quasifree ridge between 0.14 ^ ^ 0.83 (GeV/c)^ shows that there is a significant. 
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nearly constant, continuum yield extending to very large missing energies even when the 
missing momentum is relatively small. The continuum yield at large Em is predominantly 
transverse and cannot be reproduced by calculations of multiple scattering in the final state. 
A review of these data may be found in Ref. 0. One possible interpretation of these 
results is that multinucleon absorption of the virtual photon enhances the yield at large 
missing energies even when the missing momentum is relatively small. Since the present 
model does not include multinucleon absorption mechanisms, perhaps it is not surprising 
that DWA calculations underestimate the semi-inclusive cross section for ^^C(e,e'p). For 
u; = 215 MeV the ratio between experiment and calculation is approximately constant, 
but for uj ^ 650 MeV, the discrepancy increases with A. Although no calculations of 
multinucleon absorption are available for these kinematics, it seems reasonable to assume 
that such effects can become more important as either Aor u increases such that more energy 
is available to be shared among more nucleons at greater average density. Furthermore, 
Lourie et al. have suggested that kinematic focussing of the multinucleon phase space 
would exacerbate the artificial enhancement of semi-inclusive transparency measurements 
as increases. 

V. DISCUSSION 

Most previous calculations of nuclear transparency have been based upon eikonal models. 
For example, Pandharipande and Pieper [|12| studied nuclear transparency using a correlated 
Glauber model in which 

r = i|dVp,(r')PT(r') (30a) 
Pt{y') = exp |- j^^ dz" [^p„(r', r")ap„(g, p(r"))Pn(r") + ^pp(r', r'Vpp(9, p(r"))Pp(r")]} (30b) 

where Pt{t^') represents the probability that a proton struck at position r' will emerge 
without rescattering, dpN^q, pir")) represents the effective pN cross section evaluated in 
local density approximation, and 

gpNir', r") ^ gpNipo, W - r"\) 

represents the pair distribution function evaluated at central density. Although this cal- 
culation agrees well with the MIT data, it should be noted that Eq. (^) applies to Ty 
whereas the experiment measured T_i_. It is also important to recognize that despite the 
differences between the formulations of these models, the EEI model includes essentially the 
same physics. Recall that the EEI parametrization was originally based upon Brueckner- 
Hartree-Fock calculations of the effective interaction in nuclear matter which include both 
short-range correlations and Pauli blocking, but that the parameters were adjusted to im- 
prove the fit to experimental data. The correlated Glauber model includes Pauli blocking in 
its effective cross section, whereas the EEI model includes it as a density-dependent damping 
of its imaginary central interaction. Similarly, the anticorrelation between identical nucle- 
ons is represented by the pair distribution function in the correlated Glauber model or by 
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the density-dependent short-ranged repulsion in the real central EEI interaction. Pandhari- 
pande and Pieper find that Pauli blocking, effective mass, and correlation effects all play 
important roles in their transparency calculations. 

However the validity of the Glauber model is questionable at low energies where rectilin- 
ear propagation is a poor approximation, particularly for energies as low as 180 MeV used 



in the MIT experiment and the corresponding calculations of Ref. |jT2[. The Glauber ap- 
proximation also uses a zero-range approximation to the nucleon-nucleon interaction, which 
does not apply at intermediate energies even if Pauli blocking and correlation corrections are 
made. The local density approximation provides a more realistic description of the radial 
dependence of the optical potential, but includes the Glauber approximation as a limiting 
case. For proton elastic scattering the Glauber model is generally considered adequate for 
Tp ^ 800 MeV |T^,|T5|. Although integrated quantities, such as nuclear transparency, are 



probably less sensitive to the details of final state interactions, the lowest energy at which the 
Glauber approximation to nuclear transparency can be employed has not been established. 
On the other hand, the Glauber approximation is more efficient than the distorted wave 
approximation and should become sufficiently accurate for Tp ^ 800 MeV. Furthermore, 
we find that the variation among optical potentials also decreases at higher energies where 
the Glauber approximation becomes applicable. Therefore, Glauber calculations are more 
appropriate for Q"^ > 2.0 (GeV/c)^. 

Neither the optical nor the Glauber model consider in detail the distribution of absorbed 
fiux among final states of the residual system. To the extent that nucleon knockout dom- 
inates the absorption cross section at intermediate energies, most of the fiux leaving the 
elastic channel leads to multinucleon final states in the continuum which should be excluded 
from the semi-inclusive cross section used to measure nuclear transparency. Of course, it is 
much easier to exclude the multinucleon continuum theoretically than experimentally. In- 
elastic processes within final-state interactions that result in single-hole states of the residual 
nucleus would require a coupled-channels model for detailed analysis and should be included 
in the semi-inclusive cross section, but represent a rather small fraction of the absorption 
cross section. Therefore, we expect that most inelastic processes reduce transparency as 
stipulated by the optical model and the distorted wave approximation. 

Nuclear transparency for (e, e'p) at large is also usually interpreted using a Glauber 
model. As developed by Nikolaev et al. [0,^], the Glauber transparency functions for 
parallel and quasiperpendicular kinematics take the form 

= A^^pN^ jd%{l- exp [-a,o,{pN)t{h)]} (31a) 
= ]4^-p^ / {1 - exp [-aUpNm]} (31b) 



t{b) = / dznA{b,z) (32) 



where 



is the optical thickness at impact parameter b through nuclear density nA{b,z). These 
results for nuclear transparency in the Glauber model are distinguished from our definitions 
by the superscript G. The transparency for parallel kinematics, TP, is governed by the 
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total proton- nucleon cross section, atot{pN), which represents the loss of flux in the forward 
direction (parallel to the momentum transfer). The transparency for quasiperpendicular 
kinematics is larger because integration over the full angular range recaptures the part of 
the flux that corresponds to elastic proton-nucleon flnal-state interactions and hence is 
governed by the inelastic proton-nucleon cross section, ain{pN). 

Nikolaev et al. |Q demonstrated that the difference between Tjj*^ and is quite large 
at low where the inelasticity of the nucleon-nucleon interaction is relatively small. Fur- 
thermore, averaging over the acceptance for experiment NE18, they obtained an accurate 
description of the data that is intermediate between and T^^ . However, for < 1 
(GeV/c)^ the difference between 7j| and 71 is much smaller in the distorted wave approxi- 
mation than in the Glauber model (see the end of Sec. [11 A| ). Since the absorptive content 
of the optical potential for energies below the pion production threshold is determined by 
quasifree elastic nucleon-nucleon scattering, one expects T^^ ^ 7[| 7j_ for energies large 
enough to neglect Pauli blocking. On the other hand, the most important reason why the 
Glauber model calculations of Nikolaev et al. |^ produce larger transparencies than our 
optical models for ^ 1 (GeV/c)^ can be traced to differences between the deflnitions 
of the semi-inclusive cross sections used to derive and 71 that arise from differences 
between a high-energy versus a low-energy viewpoint. 

The derivation of Eq. ( ^lb| ) requires a closure sum over all flnal states of the residual 
nuclear system that contain A — 1 nucleons {e.g. Ref. [^]). Thus, the semi-inclusive cross 



section used for includes processes in which one or more secondary nucleons is ejected 
from the residual nucleus by flnal-state interactions with the primary ejectile even if the 
resultant missing energy is outside the experimental acceptance. From a high-energy view- 
point, such processes are described as incoherent elastic rescatterings of the ejectile and are 
driven by the elastic proton-nucleon cross section, aei{pN). In fact, the only flnal states 
excluded from the semi-inclusive cross section for 7^ are those which contain additional 
particles, typically one or more pions for modest u. Therefore, in the absence of Fermi mo- 
tion, 7^ for ejectile energies below the pion production threshold would be unity because 
the proton-nucleon interaction is then elastic. According to Table 1 of Ref. [^], ^ 0.97 
for ^^C or 7^ ^ 0.84 for ^osp]^ remains large even for ~ 1 (GeV/c)^. For larger ejectile 
energies decreases as the proton-nucleon interaction becomes increasing inelastic and 
approaches T^^^ for large Q^. 

By contrast, the semi-inclusive cross section used to obtain 71 within the optical model 
excludes multinucleon emission arising from flnal-state interactions. From a low-energy view- 
point, nuclear inelastic scattering which changes the internal state of the residual nucleus 
is described as a loss of flux from the elastic channel. Thus, inelastic scattering involving 
a small energy transfer between the ejectile and the residual contributes to the absorptive 
(imaginary) potential even if the flnal state remains a low-lying excitation of the {A — 1)- 
system. Since some of these low-lying flnal states fall within the missing-energy acceptance 
and are thereby included in the experimental semi-inclusive cross section, optical model 
calculations using potentials fltted to elastic nucleon-nucleus scattering probably slightly 
overestimate the loss of flux. More accurate models which account for this effect would re- 
quire coupled-channels calculations based upon optical potentials constructed for the model 
space included within the semi- inclusive cross section, but, fortunately, since the imaginary 
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part of the optical potential for intermediate energy nucleons with 100 ^ Tp < 800 MeV is 
dominated by quasifree knockout processes, this error is small and decreases as the ejectile 
energy increases. Nevertheless, it is clear that the more restrictive summation over final 
states in the optical model definition of nuclear transparency must result in 71 < T^. Fur- 
thermore, the optical model definition of nuclear transparency is more relevant to exclusive 
(e, e'p) measurements of the missing-momentum distributions for discrete final states and 
to semi-inclusive experiments of the type performed by MIT and illustrated in Fig. |^, for 
which the Glauber model would give a result, 7^ ^ 1, that is much larger than the data 
because the ejectile energy is below the nucleon-nucleon inelasticity threshold. Although 
several other assumptions employed in its derivation also fail and preclude its application 
below about 1 (GeV/c)^, the most important deficiency of this version the Glauber model for 
low is that its summation over final states is much more inclusive than the experimental 
definition of the semi-inclusive cross section. 

The Glauber model for emphasizes the inelasticity of the nucleon-nucleon interac- 
tion, whereas the optical model for 7^ emphasizes the inelasticity of the nucleon-nucleus 
interaction. Those final-state interactions which eject one or more low-energy nucleons but 
leave the missing energy within the integration range of the experiment reduce 71 but do 
not reduce T^^^. On the other hand, elastic proton-nucleon interactions which increase the 
missing energy too much reduce 7^^^ but do not reduce 7^. Therefore, the experimental 
definition of the semi-inclusive cross section for (e, e'p) reactions lies between those used for 
and 71. For relatively low the optical model is most appropriate, but will generally 
underestimate the transparency whereas the Glauber model for 71 will substantially overes- 
timate the transparency. For large where the reactive content of the optical potential is 
dominated by nucleon-nucleon inelasticity we would expect 71 to approach 7]^ from below, 
but the Glauber model is clearly much more efficient computationally than the optical model 
under these conditions. In between we would expect these two models to bracket the data. 
However, neither approach includes multinucleon absorption of the virtual photon upon a 
correlated cluster, which also might increase the experimental cross section with respect to 
direct knockout models. 

Experimentally it would be interest to observe multinucleon knockout by electron scat- 
tering and to study the kinematic dependencies of various multinucleon channels in detail. 
Although it is not possible to separate the various processes which lead to the same final 
state in a model independent fashion, the kinematic differences among them can be use- 
fully analyzed in the context of a model. For example, single-nucleon knockout from a 
deeply bound orbital leaves the residual nucleus in a highly excited state which may decay 
by particle emission and hence produce multinucleon emission without need of final-state 
interactions. Multinucleon processes of this type are probably close to isotropic in the rest 
frame of the residual nucleus and from an optical-model viewpoint should be included in the 
semi-inclusive cross section because final-state interactions are not required. On the other 
hand, one might expect secondary nucleons knocked out by interactions with the ejectile 
to appear preferentially in the forward hemisphere; such events should be included in the 
cross section for but excluded for 71. Analyses of this type require measurement of the 
angular correlation between the high-energy primary proton and one or more low-energy 
secondary nucleons. Similarly, although it is not possible to unambiguously identify events 
arising from multinucleon absorption of the virtual photon upon a correlated cluster, an 
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enhancement of the probabihty for observing two nucleons corresponding to missing mo- 
menta ±Pm with Pm ^ kp can be interpreted within a model as a signature of short-range 
correlations. Although one would expect the longitudinal/transverse character of multinu- 
cleon emission arising from the decay of a deep-hole state or from final-state interactions to 
remain consistent with the direct knockout model, multinucleon absorption of the virtual 
photon could substantially alter the structure of the response functions. In fact, there is 
some evidence at low that nucleon knockout at large missing energy is enhanced by a 
process, perhaps multinucleon absorption, that is largely transverse. However, clarification 
of these issues will require considerably more work, both experimental and theoretical. 

It may be possible to refine optical model calculations of nuclear transparency by applying 
the statistical multistep direct reaction theory of Feshbach, Kerman, and Koonin to 
evaluate the energy and angular distributions of protons which suffer final state interactions 
and thereby to estimate the fraction of the flux described by the optical model as absorption 
that actually remains within the experimental acceptances. However, it would then also be 
necessary to test those calculations against inclusive data for (p, p') . Such data is available at 



low energies, e.g. Refs. |]87| , |88[] , but is not available for Tp > 200 MeV. If data were available, 
it should also be possible to estimate the necessary corrections using a convolution procedure. 
However, these possibilities lie well beyond the scope of the present work. 



VI. CONCLUSIONS 

We have used the distorted wave approximation to evaluate nuclear transparency to 
intermediate-energy protons in semi-inclusive (e, e'p) reactions. We compared calculations 
using density-dependent effective interactions from the EEI model, which is fitted to proton 
elastic and inelastic scattering data, and the IA2 model, which is derived from a relativistic 
boson exchange model, with global optical potentials from Dirac phenomenology. For low 
energies we also considered several traditional nonrelativistic optical models. We demon- 
strated that nuclear transparency in (e, e'p) reactions is well correlated with the proton 
absorption and neutron total cross sections calculated using these models. The IA2 model 
was found to give the most accurate predictions for neutron total cross sections at Tp > 200 
MeV. The EEI model also provides accurate predictions and extends to lower energies, but 
slightly overestimates the neutron total cross sections for heavy nuclei at Tp > 300 MeV. 

For ejectile energies near 200 MeV we find that there is considerable sensitivity to the 
choice of optical model. Global optical potentials from Dirac phenomenology yield nuclear 
transparencies that are much smaller than the data, and hence are likely to overestimate 
spectroscopic factors for discrete states. Larger transparencies are obtained from the IA2 
model, but those calculations remain significantly lower than the data. Calculations using 
the EEI model predict larger transparencies than any other model considered, but still 
remain 5-10% below the data. The EEI model predicts greater transparency than other 
models for Tp < 200 MeV, but low-energy semi-inclusive data are presently lacking. At 
larger ejectile energies the sensitivity to the choice of optical model is reduced, with all 
models considered producing similar results for Tp ^ 500 MeV. Nevertheless, the calculated 
nuclear transparencies remain substantially below the NE18 data for ^ 1 (GeV/c)^, 
with the discrepancy increasing with A. Multinucleon contributions to the continuum may 
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enhance the semi-inclusive cross section for knockout and thereby artficially enhance the 
measured nuclear transparency. These contributions are likely to increase with both mass 
and energy, but quantitative estimates are not available. 

The present model is well suited to the investigation of intermediate-energy proton knock- 
out to discrete states of the residual nucleus. The dependence of attenuation factors for 
valence orbitals upon ejectile energy can then be used to investigate nuclear transparency, 
where differing radial localizations can help discriminate between interior and surface prop- 
erties of the optical potential. The asymmetry between parallel and antiparallel kinematics 
for individual orbitals can be investigated also. Recoil polarization may also provide addi- 
tional insight into the final-state interactions, and in particular may help to discriminate 
between single-nucleon and multinucleon contributions to the continuum. 

We have also examined in some detail the difference between Glauber and optical model 
calculations of nuclear transparency for quasiperpendicular kinematics. By using a much 
more inclusive summation over final states, the Glauber model emphasizes the inelasticity 
of the nucleon-nucleon interaction, whereas with a more restrictive definition of the semi- 
inclusive cross section the optical model emphasizes the role of miclcon-nucleus inelasticity. 
Therefore, the Glauber model produces larger transparency factors than the optical model, 
with the difference between the two approaches becoming quite large for < 1 (GeV/c)^. 
However, the experimental definition of the semi-inclusive cross section usually lies between 
these extremes. Although the optical model is expected to underestimate the experimental 
semi-inclusive cross section, it is much more appropriate for low than the Glauber model. 
For large where nucleon-nucleon inelasticity accounts for a much larger fraction of the 
nucleon-nucleus inelasticity, the two models should produce similar results, but the Glauber 
approximation is computationally much more efficient. 
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FIGURES 
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FIG. 1. Proton absorption cross sections for several optical models are compared with data. 
Note that for ^^C and ^°Ca the IA2 and EEI calculations sometimes almost coincide and the EEI 
calculations sometimes obscure data points. 
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FIG. 2. Neutron total cross sections for several optical models are compared with data. Note 
that for ^^C and ^°Ca the IA2 and EEI calculations sometimes almost coincide. 
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FIG. 3. Nuclear transparency data for Tp pa 180 MeV are compared with calculations of 71 
using several optical models. The data (solid points) were obtained using a 780 MeV electron 
beam, to = 215 it 20 MeV, and q ~ 605 MeV/c. Calculations were performed for selected closed 
sub-shell nuclei with mass numbers indicated by crosses. The solid line employs the EEI, dashes 
the EDADl, dots the Schwandt, and dash-dots the Comfort and Karp (CK) potential. 
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FIG. 4. Nuclear transparency data for = 1.04 (GeV/c)^ are compared with calculations of 
T± using several optical models at Tp = 650 MeV. Calculations were performed for selected closed 
sub-shell nuclei with mass numbers indicated by crosses. The solid line employs the EEI, dashes 
the EDADl, dots the IA2 optical models. 
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FIG. 5. The energy dependence of nuclear transparency is shown for ^^C (circles), ^^Ni 
(squares), and 208p|-| (diamonds) using quasipcrpcndicular kinematics for 2.1 GeV electrons. Cal- 
culations are shown for the EEl, EDADl, and 1A2 models as solid, dashed, and dotted lines, 
respectively. 



31 



